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We study the spin-fluctuation-mediated superconducting pairing gap in a weak-coupling approach 
to the Hubbard model for a two dimensional square lattice in the paramagnetic state. Performing 
a comprehensive theoretical study of the phase diagram as a function of filling, we find that the 
superconducting gap exhibits transitions from p- wave at very low electron fillings to d x 2 ^ 2 -wave 
symmetry close to half filling in agreement with previous reports. At intermediate filling levels, 
different gap symmetries appear as a consequence of the changes in the Fermi surface topology and 
the associated structure of the spin susceptibility. In particular, the vicinity of a van Hove singularity 
in the electronic structure close to the Fermi level has important consequences for the gap structure 
in favoring the otherwise sub-dominant triplet solution over the singlet d- wave solution. By solving 
the full gap equation, we find that the energetically favorable triplet solutions are chiral and break 
time reversal symmetry. Finally, we also calculate the detailed angular gap structure of the quasi¬ 
particle spectrum, and show how spin-fluctuation-mediated pairing leads to significant deviations 
from the first harmonics both in the singlet d x 2_ y 2 gap as well as the chiral triplet gap solution. 

PACS numbers: 74.72.-h,74.20.Rp,74.25.Dw 


I. INTRODUCTION 

The Hubbard model for electrons in metals is consid¬ 
ered by many to contain the essential ingredients for high- 
temperature unconventional superconductivity . 1 ® Most 
recently, it was also realized with ultracold atoms in op¬ 
tical latticesP 8 In contrast to bulk materials, the inter¬ 
action strength and e.g. the concentration of fermions 
can be changed relatively easily in optical lattice sys¬ 
tems, allowing for systematic studies of the regimes of 
doping also far away from half filling (one electron per 
site). Theoretically, there has been significant progress 
in understanding unconventional superconducting insta¬ 
bilities driven by repulsive interactions, particularly in 
the weak-coupling limit B3QH and within various numeri¬ 
cal techniques discussed in Ref. 13] and more recently in 
Refs. 1T4HT71 At present, however, the superconducting 
phase diagram of the two dimensional Hubbard model re¬ 
mains largely unknown, especially at intermediate fillings 
and at low temperatures. 

Superconductivity mediated by antiferromagnetic 
(AF) spin fluctuations near half filling was studied e.g. in 
the early work of Scalapino et a/P for a three dimensional 
paramagnetic system close to an AF instability. There 
it was found that the Coulomb repulsion between two 
electrons may give rise to a substantial Cooper-pairing 
strength due to the proximity of the AF instability. In 
this case, the superconducting gap symmetries are in¬ 
timately related to the structure of the paramagnetic 
spin susceptibility, and therefore sensitive to the geom¬ 
etry of the Fermi surface. In the work of Scalapino et 
al , 3 ^ it was found that the d x 2_ y 2 solution is favored 
close to half filling as a consequence of the spin sus¬ 


ceptibility peak at the AF wave vector Q = (7r, 7r, 7r), 
whereas the limit of very small electron filling prefers a 
spin triplet p -wave solution due to a susceptibility peak 
near q = (0,0,0). These arguments also carry over to 
two dimensions, where the corresponding wave vectors 
are ( 7 r, 7 r) and (0, 0), respectively. Since the seminal work 
of Ref. [3] spin-fluctuation pairing in the two dimensional 
Hubbard model has been studied extensively by analyti¬ 
cal and numerical approaches, and the sensitivity to the 
underlying Fermi surface gives rise to a particularly rich 
phase diagram with many emerging gap symmetries. In 
particular, it has been noted that there is support for 
triplet superconductivity in a quite large region of dop¬ 
ing in the limit of small U In Ref. [12] it was found 
that the vicinity of a van Hove singularity near the Fermi 
level causes enhancement of the triplet superconductivity 
compared to singlet solutions for sizable values of nearest- 
neighbor hopping integrals. Furthermore, it was found 
that the superconducting phase diagram in the single¬ 
band Hubbard model (with t' = 0) is quite robust against 
an inclusion of the long-range Coulomb interaction at 
least in the weak-coupling limit. 19 

In fact, the dominant interest of the community re¬ 
garding possible superconducting channels and angular 
variations of the superconducting gaps on the Fermi sur¬ 
face was often restricted to the doping near the half¬ 
filling because of the high-T c cuprates. In this regard, 
there is a consensus that the dominating gap symmetry 
of moderately hol e- or electron-doped cuprates is of the 
d x 2 _ y 2 -wave form. 20 Many experiments, however, find 
evidence that the quasi-particle gap does not always fol¬ 
low the simplest lowest order harmonic form as given by 
Ak = A[cos(^) — cos(A; y )]; the gap may exhibit its max- 
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imum value not at the antinodes, as would be the case 
for the leading harmonic d-wave form, but at a different 
location on the Fermi surface. In the electron-doped com¬ 
pounds, such a non-monotonic d-wave form was observed 
e.g. in Raman spectroscopy on NCCO 22 and in ARPES 
experiments on PLCC0. 23 Hn these cases, the position of 
the maximum gap value was related to the position of 
the so-called hot spots, i.e. segments of the Fermi sur¬ 
face which are connected by Q = (7r,7r). This points 
towards a pairing interaction mediated by AF spin fluc¬ 
tuations, in which case we expect the dominating pairing 
strength at Q, and consequently the largest gaps to be 
located at pairs of k and k' on the Fermi surface sepa¬ 
rated by Q. In the case of hole doped cuprates, ARPES 
experiments rep ort enhancement of the gap in the antin- 
odal regions. 21 This has been attributed to the presence 
of the pseudo-gap since mainly cuprates in the under¬ 
doped regime display a significant non-monoticity of the 
observed gap, and the antinodal gap is known to persist 
above T c “ However, one report on I^-^Sr^CuCU with 
hole doping of 15 %“ close to optimal doping, found a 
strong deviation from the [cos(Ay) — cos (k y )\ form. This 
indicates that the non-monotonicity may not be entirely 
caused by the pseudo-gap, and it is possible that spin- 
fluctuation pairing effects also cause gap enhancements 
near the antinodal regions of hole-doped cuprates. 

Quite generally, the interplay of spin fluctuations, su¬ 
perconducting gap and the underlying Fermi surface 
topology is an interesting characteristic of spin fluctua¬ 
tions mediated Cooper-pairing within the weak-coupling 
approach. In this regard, it is essential to study their 
evolution not only near half-filling but also for the entire 
phase diagram of the single-band Hubbard model. De¬ 
spite previous efforts in this direction 12 1 §1 18 I 25 I 26 J the S y S _ 
tematic knowledge on the weak-coupling superconduct¬ 
ing phase diagram for the single-band Hubbard model is 
still missing, especially for the intermediate doping range 
and t'/t > 0.5. 

In this paper, we investigate the spin-fluctuation- 
mediated pairing interaction within the RPA for a large 
set of possible singlet and triplet solutions to the gap 
equations within a weak-coupling approach for the en¬ 
tire doping range and strong asymmetry between elec¬ 
tron and hole doping. Concentrating on the less explored 
intermediate doping range, we show how the shape and 
topology of the Fermi surface play a decisive role for the 
final preferred gap symmetry. In particular, we find that 
the change of the Fermi surface topology, associated with 
the chemical potential crossing the van Hove singularity 
with a logarithmic divergence in the density of states, has 
a strong effect on the potential gap solutions in various 
symmetry channels, and favors a higher order triplet gap 
over the singlet d x 2_ y 2 solution even at significant elec¬ 
tron filling. We study how the different gap symmetry 
solutions evolve as a function of filling and next-nearest 
neighbor hopping integral t'/t, and map out the detailed 
gap structure arising directly from the spin-fluctuation 
pairing mechanism. This includes deviations from the 


[cos (Ay) — cos(fc y )] form of the superconducting gap in 
the singlet channel close to half filling, as well as the 
form of the higher order triplet gap. In the triplet chan¬ 
nel we find that the preferred solutions are those that 
break time reversal symmetry. 


II. MODEL AND METHOD 

We consider the Hubbard model for a two-dimensional 
square lattice 


H — ^^ C kcr C ka + 2 N ^ C k / (7 C -k / +qq C -k+q<r c k(r; 

k a k,k',q cr 

( 1 ) 

where £ k = — 2t [cos (Ay)+cos (k y )] —At' cos(Ay) cos {k y ) — n 
with t being the hopping integral to nearest neighbors, 
and t' < 0 the hopping integral between next-nearest 
neighbors. In the following we set t = 1 and restrict 
ourselves to the case of negative values of t'. 

A spin-fluctuation-mediated interaction can combine 
two electrons of opposite spin or the same spin into a 
Cooper pair. The pairing interaction is derived from 
higher order diagrams of the repulsive Coulomb inter¬ 
action [/EC 27 i n the case of opposite electron spins, the 
diagrams consist of an even number of bubbles as well 
as ladder diagrams, which correspond to spin preserv¬ 
ing or spin flip interactions, respectively. For same spin 
electrons, the interaction is derived from an odd number 
of bubble diagrams, and in this case only spin preserving 
interactions are allowed. Specifically, the interactions are 
given bj0 
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with the spin and charge susceptibilities given by 
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Here U = U/z is a renormalized Coulomb interaction. In 
the main part of the paper we use z = 1 corresponding 
to the usual RPA, but when solving the full gap equa¬ 
tion we use a renormalization (z = 2,3) to achieve a 
larger pairing strength for numerical convergence. Equa¬ 
tions © and Q provide a measure of the interaction 
strength, and we neglect the energy dependence of the 
interactions. The bare susceptibility in the paramagnetic 
phase is given by the Lindhard function 


xo(q,w) = 


w + a - a + q +vi ’ 


(6) 
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which is evaluated a zero energy (uj = 0). The gap equa¬ 
tion arises from a standard mean-field decoupling of the 
interaction Hamiltonian. In the singlet ( s ) and triplet ( t ) 


channel it takes the form 

1 A s E 

A s / f — 1 Vr s/ * k ' tanhl 

' E * \ 


A k - 2N ^ o F s/t ' 

k' Z£, k' 

v2 k B TJ’ 

(7) 

with 

E± = \4k + |Aifl 2 - 


(8) 


In the calculation of the superconducting gap, the po¬ 
tential forms stated in Eqs. © and © must be sym¬ 
metrized or antisymmetrized with respect to momentum 
in the singlet and triplet channel, respectively. In case 
of singlet pairing, the interaction potential Tk k 7 is given 
by the opposite spin vertex, as stated in Eq. Q, and we 
have 


ps _ popp.sp , F opp.sp 

1 k,k' — 1 k,k' 1 —k,k' * 


(9) 


For triplet pairing, the same effective interaction is ob¬ 
tained irrespective of whether the same spin or opposite 
spin interaction vertex is used 


pt _ popp.sp _ popp.sp _ psame sp _ psame sp (10) 

k,k^ k,k^ k 5 k^ k 5 k^ \ z 

This is as expected in the paramagnetic phase. Note 
that the potential entering Eq. (|t|) appears in the singlet 
(even in k) and triplet (odd in k) form explicitly. This 
symmetry directly carries over to the gap, ensuring that 
A k = Ai k and A k = — A^_ k . When solving the gap 
equation as stated in Eq. 0 , we invoke a small energy 
cut-off around the Fermi surface, e c , and allow for Cooper 
pair formation of all electronic states with £k £ [—e c , e c ]. 

In order to capture the filling dependence of the most 
prominent gap candidates, we project the pairing poten¬ 
tial onto the leading order gap harmonics given by 


s* = cos (k x ) + cos (fc y ), (11a) 

d x 2_ y 2 = cos (k x ) — cos (fc y ), (Hb) 

d xy = sin(fcaj) sin(fcy), (11c) 

g = [cos(^) — cos(fcy)] sin (k x ) sin (k y ), (lid) 

pm sm(k x ), (lie) 

p' = [cos(^) — cos (k y )] sin (k x ). (Ilf) 


Note that in the literature on Sr 2 Ru 04 the p' solution is 
sometimes called / x 2 _ y 2 - wave The two triplet solutions 
p and p' belong to the same two-dimensional E u group, 
and are visualized in Fig. [l] We follow the procedure 
of Scalapino et al. 3 and calculate the projection of the 
interaction vertex onto the basis functions 
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FIG. 1. (Color online) Illustration of the triplet gaps p-wave 
(a) and p -wave (b). The latter is favored in the doping 
regime where the spin susceptibility has a peak or plateau 
at Q = (7r,7r). The magnetic zone boundary is indicated by 
the dashed red line. Nodal lines in the two cases are shown 
by full black lines. 


with g a being one of the functions stated in Eqs. (11a)- 
( |llf| ). This procedure does not include possible higher 
order solutions. To determine these, we also solve the 
linearized gap equation in the singlet and triplet channel 


h?S, 


dk' s /t 

i i J- 1 
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2(2tt)* j F s |Uk' 
by diagonalization of the matrix 
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Here k and k' are located on the Fermi surface and l k is 
the length of the Fermi surface segment associated with 
the point k. By this procedure we identify the leading in¬ 
stability as a function of electron filling and next-nearest 
neighbor hopping constant, t'. We characterize the lead¬ 
ing singlet solution according to its transformation prop¬ 
erties into one of the four singlet representations A\ g : 8*, 
B\g : d x 2 _ y 2 , A^ g : p, B 2 g : d xy . Note that the square 
lattice has one class of triplet solution (E u ) of which the 
p -wave is the lowest harmonic. However, in general we 
find the leading triplet solution to be higher order, as 
discussed in detail below. 

The linearized gap equation does not allow for a deter¬ 
mination of complex gap solutions, which are time rever¬ 
sal symmetry broken (TRSB) solutions. TRSB solutions 
can lead to a removal of gap nodes from the Fermi sur¬ 
face whereby there is a gain in condensation energy. In 
the triplet channel, where solutions are doubly degen¬ 
erate, TRSB solutions might be favored due to this ef¬ 
fect. Therefore, we also address the non-linearized gap 
equation as stated in Eq. 0 and show that the full self- 
consistent calculation finds that the triplet solutions are 
TRSB. Finally the solution of the full gap equation is 
used to obtain the angular structure of the resulting gap 
which may exhibit significant changes from the standard 
lowest harmonic due to the details of the momentum 
structure of the spin susceptibility. 
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FIG. 2. (Color online) (a-j) Spin susceptibility in the RPA approximation and (k-t) Fermi surfaces for fillings (n) = 
0.05,0.35,0.50,0.65,0.80 and (n) = 1.05,1.2,1.35,0.55,0.90 with t' = -0.35 and U - 1.75. The AF zone boundary is shown 
by the dashed red line. Note the different colorbar range for each susceptibility plot, (u) Phase diagram for A as stated in 
Eq. (12) for a band with t' = —0.35 and U — 1.75. The temperature is ksT = 0.015. Projection onto s*-wave (cos k x + cos k y ) 
gives a negative value of A at all fillings. The left inset shows a zoom for low filling values (n) — 0 — 0.2. The inset to the right 
shows the density of states at the Fermi level, p(0), as a function of filling. Note that at the filling for which p(0) is maximal 
the triplet p' solution is the leading instability. This we refer to as the van Hove critical density. For t' — —0.35 the van Hove 
critical density is (n) v #=0.66. 


III. RESULTS 


A. Doping dependence of the lowest harmonic 
solutions 


First we consider the Fermi surfaces and spin suscep¬ 
tibilities throughout the entire doping range for a next- 
nearest hopping integral t' = —0.35 relevant for cuprates. 
As opposed to the early work in Ref. [3] which reported 
the t f — 0 case, we pay special attention to the electron- 
hole asymmetric case when the next-nearest neighbor 
hopping t' is non-zero. In Fig. [2jk-t) the Fermi surfaces 
at different electron fillings are shown. Note the transi¬ 
tion of the Fermi surface between Fig. [2p) and Fig.^o). 
In the latter case, the Fermi surface has ” split up” at 
the antinodal positions (7r,0) and (0,7r). This splitting 
occurs when the van Hove singularities at (±7r, 0) and 
(0, ±7r) cross the Fermi level for a chemical potential of 
/i = At'. For t' = —0.35 this happens at the van Hove 
critical density, ( n) V H = 0.66. 

The spin susceptibilities at the corresponding doping 
levels are shown in Fig. [2| a-j). At very large dopings the 


susceptibility exhibits a broad peak around q = (0,0) as 
seen in Figs. |2^a) and Hi) which develops into peaks at 
(dz7r, 0) and (0, ±7r) as the doping is decreased as seen 
for (n) = 0.35 in Fig.[2|b). At intermediate doping levels 
the peaks at (±7r, 0), (0, ±7r) move inwards, and develop 
into the well-known quartet of incommensurate peaks at 
( 7 r dz S,7r) and (n, it db J) as shown in Fig. [2jc), with 5 
decreasing as the system gets closer to half filling. A 
special feature is observed for fillings close to the van 
Hove critical density as seen from Fig. [2|d) where a q = 
(0, 0) peak develops as a direct consequence of the large 
density of states at the Fermi level. Close to half filling, a 
clear peak around Q = ( 77 , 7 r) develops as shown Fig.[2je- 
g). Finally, at very large fillings the spin susceptibility 
becomes almost featureless as seen from Fig. Hh-j). 

In order to map out the fillin g de p ende nce of the 
gap symmetries as defined in Eqs. (lla)-( pTf| we project 
the pairing potential onto these symmetries as stated in 
Eq. (12). In the following, we do not show the s*-wave 
results since these are highly suppressed at all doping lev¬ 
els. We have chosen the Coulomb interaction U = 1.75. 
While this choice leads to small values of A, it allows us 
to avoid the instability to long range magnetic order over 
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FIG. 3. (Color online) Phase diagram of the leading super¬ 
conducting instability as a function of electron filling (n) and 
next-nearest neighbor hopping constant \t'\. Throughout, the 
value of U is adjusted such that the leading eigenvalue is 
A = 0.1. For values of U <3, the eigenvalue is shown in 
full color. For large U values in the range 3 — 8 the eigenvalue 
is shown with brighter color in order to indicate this. We in¬ 
clude only fillings (n) <1.5 since the instability is negligible 
for higher electron dopings for U <8. The three green stars 
mark the positions for which we plot the singlet solutions in 
Figs. [8| and [To] and the two black filled circles mark the posi¬ 
tion for which we solve the full (non-linearized) gap equation 
in the triplet channel, Fig. [7] 


the entire phase diagram. At very large hole and elec¬ 
tron dopings, the peaks in the spin susceptibility are weak 
in intensity and, as a direct consequence, the projected 
pairing strengths are relatively small compared to half 
filling. This doping regime has a leading triplet p -wave 
solution, which for the hole doped case is visible from the 
left inset of Fig. [2| This is in agreement with the result 
of the t' = 0 bandP In the filling regime (n) ~ 0.2 — 0.4, 
the g-wave solution is favored, whereas none of the low¬ 
est harmonic solutions are supported in the filling regime 
around (n) = 0.45. It is important to note that the pro¬ 
jection method does not take into account possible higher 
order solutions. These might dominate in some regions. 
For instance, the absence of a positive A for (n) ~ 0.45 in¬ 
dicates that the leading solution in this region is a higher 
order solution. For the more general solution we deter¬ 
mine the leading solution to be a higher order 5*-wave in 
this case, as we shall see below. 

Close to half filling the d x 2_ y 2 solution is clearly dom¬ 
inant. When the system is hole doped away from half 
filling, the d x 2_ y 2 solution becomes increasingly strong 
as the the van Hove critical density, ( u) v h = 0.66, is ap¬ 
proached. However, an abrupt change occurs very close 
to the van Hove critical density, where the d x 2_ y 2 solu¬ 
tion becomes unstable. This is manifested by a sharp 
dip of the dashed cyan curve in Fig. [2p). The d x 2_ y 2 




FIG. 4. (Color online) (a) Spin susceptibility Re xo(q, u = 0) 
at the wave vector q = (0, 0) as a function of filling, (n) and 
next-nearest hopping constant t'. The full black line and the 
black dashed-dotted line show where the van Hove singularity 
crosses the Fermi level at the positions (d=7r, 0)/(0, d=7r) and 
(k, =b k), with k = d= cos -1 (^q)? respectively. The dashed line 
shows where a hole pocket is removed from the Fermi surface. 

(b) Spin susceptibility weight at the wave vector Q = (7r,7r). 

(c) The wave vector q, for which the bare spin susceptibility 
achieve its maximum value (color indicates q according to 
inset), plotted as a function of filling (n) and next-nearest 
hopping constant t f . 


solution becomes unfavorable due to the development of 
a peak at q = (0,0) in the susceptibility. This gives 
rise to a large repulsive interaction between neighboring 
momenta k and k' in the singlet channel, and causes a 
suppression of all singlet solutions, in this case the d x 2_ y 2 
solution. Such a suppression of singlet superconductivity 
due to the q = (0,0) peak in the spin susceptibility was 
originally discussed by Berk and Schrieffer 27 in their pi¬ 
oneering work on spin-fluctuation mediated pairing. As 
clearly visible in Fig. [2|u) the sharp dip of the singlet 
solution around the van Hove critical density is accom¬ 
panied by an increase in the triplet solution p' shown by 
the red line. Thus, the q = (0, 0) susceptibility peak not 
only suppresses the singlet solution, but actually sup¬ 
ports the development of a triplet gap because it gives 
rise to an effective attraction for neighboring k and k' at 
the Fermi surface. In the absence of additional structures 
in the spin susceptibility the p -wave solution is favored, 
as in the case of very small filling shown in the inset of 
Fig. |2|u). However, if the spin susceptibility shows addi¬ 
tional peaks as in the case of fillings close to the van Hove 
critical density, Fig. |2](d), higher order triplet solutions 
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will be favored. The detailed structure of the potential 


turns out to favor a six node p'-wave gap, see Eq. ( Ilf) 
and Fig. 0b). We will return to this is section [III C| 


B. Phase diagram as a function of doping and t' 


In order to further investigate the generality of the 
results obtained by projection in Fig. [2] we turn to a 
solution of the (non-projected) linearized gap equation, 
Eq. (13), for different next-nearest neighbor hopping 


strengths \t'\ in the range 0 — 1.5. In contrast to the pro¬ 
jection method, this general procedure determines the 
leading solution amongst all possible higher order har¬ 
monics. We classify the solution according to its trans¬ 
formation properties into one of the four singlet solutions 
s*, d x 2 _ y 2 , d xy , g or triplet. In the two subsequent sec¬ 
tions, we will discuss the triplet states in more detail. 
Note that in general solutions may correspond to higher 
order harmonics and have additional nodes co mpared to 
the leading harmonics given in Eqs. ( lla|lle ). For each 
point in the phase diagram ((n), t') we solve Eq. ( fTU] ) and 
adjust the value of U such that the leading eigenvalue is 
A = 0.1. The result of this procedure is shown in Fig. [3j 
where the leading superconducting instability for fillings 
and next-nearest neighbor hopping strengths in the range 
0 — 1.5 is shown. This procedure is justified as follows. 
Since we cover a large range of different Fermi surface 
structures, a fixed value of the Coulomb interaction of 
e.g. U — 2 will cause a break-down of the paramag¬ 
netic RPA formalism due to the instability to long range 
magnetic order. By allowing for a variation in U this is 
avoided, and at the same time we discuss only instabili¬ 
ties with a non-negligible critical temperature. Note that 
this approach is different than a previous report P 8 where 
the genuinely weak-coupling approach U 0 limit was 
taken. Nevertheless, our results show some qualitatively 
agreement with Ref. [18j a g-wave and a (small) d. 


xy 


re¬ 


gion appear in the filling regime of (n) = 0.25 — 0.4 for 
\t'\ < 0.5 and at fillings around (n) = 0.7 a d x 2_ y 2 region 



FIG. 5. Fermi surfaces for t' m —0.8 at three fillings, (a) 
At small fillings, (n) =0.3, the Fermi surface consists of four 
electron pockets centered at (d=7r, 0) and (0, d=7r). (b) For a 
filling of (n) =0.9, which is above the van Hove critical den¬ 
sity, a hole pocket is centered at (0, 0) and four hole pockets 
are centered at (d=7r, d=7r) and (d=7r, =F7r). (c) At a larger filling, 
(n) = 1.14, the central hole pocket is removed from the Fermi 
surface. 


dominates at all \t'\ < 0.4, whereas an s-wave domain 
takes over for 0.4 < \t'\ < 0.5 in this filling regime. In 
our approach, the region of triplet solutions has substan¬ 
tially shrunk compared to Ref. QS with different higher 
order singlet solutions taking over in the regime of small 
t' close to a filling of (n) =0.5. In the case of large elec¬ 
tron dopings, where the spin susceptibility shows only 
very weak structure, U must be very large to obtain a 
leading instability of A = 0.1 and we omit this regime in 
Fig-! Further, regions for which U > 3 the colors are 
less saturated to indicate this. 

The structure of the spin susceptibility plays a decisive 
role for the gap symmetry, and this is very tightly con¬ 
nected to the geometry of the Fermi surface. Two differ¬ 
ent regimes of t' give rise to very different Fermi surface 
geometries. For \t'\ < 0.5 the Fermi surface evolution 
as a function of doping is similar to the case t' = —0.35 
shown in Fig. [2|k-t). In this case, a transition of the 
Fermi surface occurs when the van Hove singularities at 
(±7r, 0) and (0, ±7r) cross the Fermi level. This happens 
when fi = At'. For \t'\ > 0.5 the transition occurs at 
g = 1/t' and invokes van Hove singularities at diagonal 
positions k y = ±k x = db cos_1 ( 2 fF|)' The Fermi surface 
evolution with doping for t' = —0.8 is depicted in Fig. [5] 
For the special case of \t'\ =0.5 the van Hove singularity 
resides at the bottom of the energy band. 

The correlation between the van Hove singularity at 
the Fermi level and the appearance of spin susceptibility 
weight at the wave vector q = (0, 0) is shown in Fig.[4^a). 
Here the value of the spin susceptibility at q = (0, 0) is 
shown as a function of (n) and t'. We also show where the 
van Hove singularity crosses the Fermi level, indicated by 
the full black line in the case of \t'\ < 0.5) and the dashed- 
dotted line for the \t'\ > 0.5 case. The black dashed line 
shows where a hole pocket is removed from the Fermi 
surface. In Fig. J^b) the spin susceptibility weight at 
the wave vector Q = (7r,7r) is depicted. As expected, 
the signatures of a strong peak at Q is clearly visible 
near half filling and t! = 0. The weight at Q expands 
in a doping region around half filling for \t'\ < 0.5. In 
addition, it shows a clear correlation with the van Hove 
critical density shown by the full black line. This explains 
why the d x 2_ y 2 solution is increasingly favored upon hole 
doping away from half filling for t' = —0.35, as observed 
in Fig. 0u). 

By a comparison of the red regions in Fig. [3] and 
Fig- 0a), we observe that the two triplet solution 
branches that expand from the low filling regime in Fig. [3] 
are explained by the weight at q = (0,0) in the spin sus¬ 
ceptibility, which is correlated with the van Hove critical 
densities. This is in agreement with the expectations 
that a q = (0,0) peak in the pairing interaction favors 
triplet superconductivity. The reason why we see a shift 
from triplet to singlet ( d x 2_ y 2 ) superconductivity at the 
lower red branch around (n) ^ 0.6 is because the Q peak 
becomes dominant in this regime. This transition from 
triplet to d x 2_ y 2 superconductivity is most clearly visible 
in Fig. |4jc), where we capture the dominating structure 
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FIG. 6. (Color online) Singlet (dashed lines) and triplet (full lines) pairing interaction in the form stated in Eq. 0 . 
for filling (a) (n) = 0.80 and (b) (n) = 0.65. Two different positions of the momentum k has been chosen, as shown by the 
colored crosses in the Fermi surface inset. The pairing between k and all k' in region (1), (2), (3) and (4) of the Fermi surface 
are shown in the main panels with the position of k' being parametrized by the angle 0 measured with respect to the x axis. 
The next-nearest neighbor hopping is t' — —0.35 and the Coulomb interaction is U — 1.75. 


of the spin susceptibility at every position ((n),£') of the 
phase diagram, by plotting the wave vector q for which 
the bare spin susceptibility achieve its maximum value. 
Note how the structure of the spin susceptibility trans¬ 
forms from being dominated by a q ~ (0, 0) peak in the 
regime of the triplet branch, to being dominated by the 
Q peak. Also, from the large magenta region in Fig.[4^c) 
we see why triplet superconductivity governs the phase 
diagram in an extended regime around \t'\ = 0.5 for small 
fillings. As expected, a Q peak dominates in the region 
around half filling. 

In other regions of the phase diagram where we ob¬ 
serve a different type of singlet superconductivity, other q 
structures of the susceptibility become dominant, as seen 
in Fig. |4jc). Note, however, that subdominant features 
in the susceptibility which might influence the gap equa¬ 
tion, are not visible from this figure. The large region of 
d xy superconductivity, which occurs at small to moderate 
hole dopings and small electron dopings for \t'\ > 0.5 is 
correlated with spin susceptibility peaks near ( 7 r, 0) and 
(0,7r). We return to the different manifestations of the 
d xy for small and large \t'\ in section 


IIIF 


C. Triplet gap at the van Hove critical density 

In the section above we saw how the suppression of 
singlet superconductivity and concurrently, the develop¬ 
ment of a triplet gap is intimately related to a q = (0, 0) 
peak in the susceptibility which occurs at the van Hove 
critical density. Now we turn to a more detailed investi¬ 
gation of the structure of the pairing potential and the 
consequences for the favored gap symmetries. In our 
model two contributions are important for the pairing 
structure: 1) the hot spot effect, which for most filling 


levels is more accurately described as a plateau around Q 
rather than a sharp peak at Q as seen e.g. from Fig. [2^f ), 
2) the van Hove effect with pairing contributions arising 
due to the appearance of a quartet of peaks at the di¬ 
agonal corners at q v h = (d, d) where d 0 as the van 
Hove singularity crosses the Fermi level (in the case of 
\t'\ < 0.5). The peaks at q v h are visible in Fig. |2|e) and 
appear as a purple region in Fig. [4jc) when hole doping 
is increased towards the van Hove critical doping, shown 
by the full black line in Fig. j4^c). To visualize the hot 
spot effect and the van Hove effect explicitly, we plot the 
pairing potential for the band with t' = —0.35 and fill¬ 
ings (n) = 0.80 and ( n ) = 0.65 in Fig. [6ja) and|6|b), 
respectively. In the case of (n) = 0.80, the largest pair¬ 
ing potential is found due to the hot spot effect, as seen 
by the red curve close to the angles 0 = | and 0 = ^ 
in Fig. [(^a). A smaller signature due to the van Hove 
effect is seen in the red curve close to 6 = 0 and 0 = tt. 
As the hole doping is increased, the van Hove effect be¬ 
comes more pronounced, and very close to the van Hove 
critical density at ( n) V H = 0.66, we observe sharp peaks 
at 6 = 0 and 0 = n in the red dashed curve in Fig. |6|b). 
These peaks are responsible for suppression of the singlet 
solution. At the same time the van Hove effect gives rise 
to attractive potentials in the triplet channel, as seen in 
from the full blue line at 0 = j and the full red line at 
(9 = 0. 


At very small fillings, the triplet solution has the sim¬ 
ple p-wave form as stated in Eq. (lie), but away from the 
small filling regime, the triplet solution is represented by 
higher order harmoni cs. O ne of these solutions is the p' 
solution given in Eq. ( |llf| ). The main difference between 
the p-wave and the p'-wave is that the latter gap has the 
same sign for sgn[Ak + q] = sgn[Ak], whereas the p-wave 
solution obeys sgn[Ak+Q] = — sgn[Ak]. These proper- 
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FIG. 7. (Color online) (a) Solution to the full gap equation 
Eq. Q with triplet potential and the bare Coulomb interac¬ 
tion U = 7.4 renormalized to U — 3.7. The filling level is 
(n) = 0.03. Real (full red line) and imaginary (black dashed- 
dotted line) part of the chiral p -wave gap as a function of angle 
6 defined in the Fermi surface inset, (b) Solution to Eq. 0 
with triplet potential and bare Coulomb interaction U = 5.7 
renormalized to U — 1.9. The filling level is (n) — 0.65. Real 
(full red line) and imaginary (black dashed-dotted line) part 
of the chiral higher order triplet gap as a function of angle 0 
defined in the Fermi surface inset. The absolute value of the 
gap |A(0)| is shown by the blue dash dotted line. In both 
cases, the next-nearest neighbor hopping is t' = —0.35, tem¬ 
perature is ksT — 0.01 and energy cut-off is set to e c = 0.012. 
The size of e c does not change the gap structure qualitatively. 


ties are illustrated in Fig. [T] When the susceptibility has 
a peak or a plateau at Q = ( 7 r, 7 r), it is favorable for 
the triplet gap to display the same sign at k and k' dis¬ 
placed by Q, since the triplet pairing potential contains 
the term — ^rx(k — k'). Therefore the p'-wave gap will 
be favored in the filling regime where the susceptibility 
shows a q = (0,0) peak as well as a peak or plateau 
structure around Q = (7r, 7 r). In the more general case of 
a quartet peak structure around Q as seen in Fig. [2|d) 
the triplet gap will resemble the structure of p', but with 
the nodes slightly displaced, which we will discuss in the 


next section. 

In Ref. |32| the changes in gap symmetry as a function 
of electron- and hole doping within spin fluctuation me¬ 
diated pairing to second order in U = 6t were discussed. 
The pairing potential was H(k, k') = U + U 2 x{ k + k'), 
and the next-nearest hopping constant fixed at t'/t = 
—0.276. The d x 2_ y 2 solution dominates at all moderate 
doping levels, also at the van Hove critical density, in 
agreement with our findings. However, in contrast to our 
results, Ref. [32] reports a regime of triplet superconduc¬ 
tivity at smaller fillings. As mentioned above, Ref. 18 
also reports an extended region of triplet superconduc¬ 
tivity which is unrelated to the van Hove critical density. 
We suspect this discrepancy to arise from the details of 
procedure; in Refs. Ql| and 32] gap solutions was trun¬ 
cated to the first 15 harmonics, whereas in our case we 
do not invoke any restrictions on the gap functions. 

In a recent work by Deng et al. the emergence of 
pairing for the paramagnetic liquid was also addressed 
in a numerical study of the two-dimensional Hubbard 
model with i! m 0. They reported a transition from p- 
wave superconductivity at small fillings through a d xy 
gap at intermediate filling levels to a d x 2_ y 2 symmetry 
close to half filling. For small values of 7/, a higher order 
triplet gap with six nodes was also found for fillings at 
(n) ~ 0.55. This triplet solution is thus unrelated to the 
van Hove critical density, and we do not find a similar 
solution at t' = 0 in our case. Their findings, however, 
agrees with the results in Ref. [18] which report the U 0 
case. The disagreement might arise from the difference 
in strength of U . 


D. Time reversal broken triplet gap solutions 

The fact that all triplet solutions found in the lin¬ 
earized approach are two-fold degenerate suggests that 
a TRSB solution might be favored. Therefore we turn 
to the full gap equation as given in Eq. 0 in the triplet 
channel, which we solve at ksT = 0.01. We consider two 
special fillings of (n) = 0.03 and (n) = 0.65 where in 
both cases t' = —0.35. At very low filling, the suscep¬ 
tibility exhibits only a weak structure around q = (0,0) 
as seen in Fig. [2|a). In this case the p -wave triplet so¬ 
lution is favored and, as shown in Fig. J7^a), the full so¬ 
lution is the nodeless TRSB gap of the form p x db ip y . 
In this filling regime, strong Coulomb interactions are 
required to achieve a superconducting instability due to 
the weak structure of the spin susceptibility. We use a 
bare Coulomb interaction of U = 7.4 which is renormal¬ 
ized in the RPA expressions to U = 3.7. 

At higher fillings, the spin susceptibility acquires more 
structure and supports a superconducting gap with a 
bare Coulomb interaction U = 5.7 renormalized to U — 
1.9. The preferred solution in this case is also a TRSB 
solution. If the TRSB solution had been of the form 
p' x dz ip y only the nodes along the zone axes k x or k y 
would be lifted with the four nodes along the zone di- 
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FIG. 8. (Color online) (a) d x 2_^2 solution to the linearized 
gap equation as stated in Eq. (fl3j) for a band with t' = —0.35 
and U = 1.75 at a filling of (n/ = 0.8. (b) d x 2_ v 2 solution to 
the linearized gap equation as stated in Eq. (13) for a band 
with t' = —0.8 and U = 1.75 at a filling of (n) = 0.5. The 
magnetic zone boundary is indicated by the red dashed line. 



0[tt] 


agonals preserved. In this case, since only two of the 
six nodes are lifted, the gain in condensation energy of 
the TRSB solution compared to one of the solutions, p' x 
or p f y , would be limited. However, in the present case 
of (n) = 0.65 where the spin susceptibility displays a 
quartet of peaks around Q, the simple form p' x =t ip' is 
replaced by a more complicated gap solution displayed 
in Fig.[7jb), which indeed provides a fully gapped TRSB 
solution. From the angular gap dependence shown in 
Fig- 0b), strong effects due to the susceptibility struc¬ 
ture are clearly visible from the absolute value of the gap, 

I Atrip (0)| as shown by the full blue line. The maximum 
of the gap achieved close to 0 = 0 and | is due to the 
van Hove effect, and the peaked feature at 6 ^ i r/8 is 
related to the quartet of peaks around the Q vector in 
the susceptibility as shown in Fig. [2jd) . This underlines 
again the strong connection between the Fermi surface 
structure, spin susceptibility, and the detailed angular 
dependence of the superconducting gap. Lastly, we note 
that the gap minimum of this higher order triplet solu¬ 
tion is achieved at 0 = J, which is at the position of the 
nodal lines of the d x 2_ y 2 structure. 


E. The d x 2_ y 2 solution 

The d x 2_ y 2 solution of the one-band Hubbard model is 
commonly discussed in the region with \t'\ < 0.5 around 
half filling, due to the relevance for cuprates. From the 
phase diagram in Fig. [3j we observe that there is in fact 
another large region of the phase diagram for which a 
d x 2_ y 2 solution dominates corresponding to \t'\ > 0.5 
where the Fermi surface topology is quite different. In 
order to show a solution in both regimes, we plot in Fig. [8] 
the d x 2_ y 2 solution in the case of t' = —0.35, (n) = 0.8 
and t' = —0.8, (n) = 0.5. Note that even though the 
Fermi surface of the t' = —0.8 band has no Fermi surface 
weight along the zone diagonals, which are the nodal lines 
of the simple harmonic d- wave, it displays more nodes at 
the Fermi surface than the solution for the t' = —0.35 


FIG. 9. (Color online) Solution to the full gap equation 
Eq. § with singlet pairing interaction for t' — —0.35 at fill¬ 
ing levels (n) =0.85 (full red line) and (n) m 0.75 (dashed 
blue line) at temperature ksT — 0.01. The renormalized 
Coulomb interaction is U = 2,1.88, respectively and the 
Coulomb renormalization is z = 3. The energy cut-off is set to 
e c = 0.015, and the gap value at the antinodes is A^> =0 = 0.01 
in both cases. The angle 0 is defined in the Fermi surface in¬ 
set. The vertical lines show the position of the hot spots for 
the (n) = 0.85 (full red line) and (n) = 0.75 (dashed blue 
line). Note that there is a clear correlation between the posi¬ 
tion of the hot spot and enhancement of the superconducting 
gap away from a simple harmonic gap function in the case 
of (n) = 0.85. However, at (n) = 0.75, which is closer to 
the van Hove critical density, the higher harmonic angular 
dependence of the gap function is not directly related to the 
position of the hot spot. 


band, namely eight nodes instead of four. 

A closer inspection of the d x 2_ y 2 solution for the band 
with t' = —0.35 indicates deviations from the simple form 
of A [cos k x — cosky] due to the presence of higher har¬ 
monics, i.e. longer range superconducting pairing inter¬ 
action. The maximum gap value is not achieved at the 
antinodal points, but shifted towards the nodal direction. 
This effect has been discussed previously as a signature 
of spin-fluctuation-mediated pairing!j n p rev i 0 us 
work special attention was drawn to the hot spot effect 
in which the gap maximum occurs at the k position of 
the hot spot. In Fig. [9] we show the singlet gap as calcu¬ 
lated by the full (non-linearized) gap equation, Eq. (|7|). 
In the figure we also show the position of the hot spot, 
i.e. the angle at which the Fermi surface intersects the 
magnetic zone boundary. This is shown by the verti¬ 
cal lines in Fig. [9j It is seen that the non-monotonicity 
of the gap is not directly related to the hot spot effect, 
since the hot spot position moves towards 0 = 0, but the 
strong gap enhancement moves closer to the nodal direc¬ 
tion upon increased hole doping. Close to the van Hove 
critical density, the higher harmonic content of the gap 
function becomes more pronounced. This tendency was 
also pointed out in Ref. [32j Note that proximity of the 
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FIG. 10. (Color online) Leading solutions to the linearized 
gap equation as stated in Eq. ( |13| ) for a band with t' — —0.35 
and U — 1.75. (a) At filling of (n) = 0.25 the singlet g- 

wave solution dominates, (b) At intermediate filling level of 
(n) =0.50 the singlet higher order s-wave solution with eight 
nodes is the leading instability. The magnetic zone boundary 

is indicated by the red dashed line. 
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FIG. 11. (Color online) (a) dxy solution to the linearized gap 
equation as stated in Eq. (13) for a band with t' = —0.35 
and U = 1.5 at a filling of (n) = 0.55. (b) d xy solution to 
the linearized gap equation as stated in Eq. |l3| for a band 
with t' = —1.2 and U = 1.5 at a filling of (rJy== 0.65. The 
magnetic zone boundary is indicated by the red dashed line. 


regime of |P| >0.5 the d xy solution is largely dominating 
in a large region close to half filling. In Fig. [TT] w e show 
the d xy solution in the case of \t'\ < 0.5 (Fig. |ll[a)) and 
\t'\ > 0.5 (Fig. [llpo)). Due to the difference in Fermi 
surface topology, the appearance of the d xy solution is 
quite different in the two cases. 

IV. CONCLUSIONS 

In this paper we have studied the superconducting 
gap structures in a single band Hubbard model within 
spin-fluctuation-mediated Cooper-pairing scenario in the 
weak-coupling paramagnetic limit for an extended re¬ 
gion of phase space. It complements our earlier study 
of pai ring i n the spin density wave phase in the same 
modelf 34 35 In contrast to previous studies of the param¬ 
agnetic phase, our main emphasis was to study the gap 
structure for a large range of next-nearest neighbor hop¬ 
ping integrals, £', and doping levels away from half-filling, 
which could be potentially relevant for future systems in¬ 
cluding new classes of unconventional superconductors as 
well as optical lattices loaded with interacting fermions. 
We discussed the details of the gap structure and related 
this directly to the spin susceptibility at all filling lev¬ 
els. Furthermore, we also focused on the role of a van 
Hove singularity in close proximity to the Fermi level 
for the transition between various Cooper-pairing chan¬ 
nels. This has drastic effects on the gap symmetry since 
it strongly suppresses singlet superconductivity and leads 
to the emergence of a nodeless TRSB triplet gap solution. 
This is a direct consequence of the additional q = (0,0) 
peak structure in the spin susceptibility which reflects 
the presence of a van Hove singularity at or very near to 
the Fermi level. 


van Hove singularity leads to an increase in the number of 
states participating in the formation of Cooper pairs, but 
the corresponding formation of additional nesting peaks 
in the susceptibility at small = (S,S), see Fig. [2je), 
in fact work against the d x 2_ y 2 solution, since in the sin¬ 
glet channel this will favor nodes along the zone axes. 


F. Other singlet solutions 


From the phase diagram in Fig. [3] we observe that the 
regime of \t'\ < 0.5 and intermediate hole doping levels, 
have two robust regions of singlet superconductivity be¬ 
sides the d x 2_ y 2 solution, namely a g-wave and s-wave 
region. We show the solutions for t f = —0.35 at fillings 
(n) =0.25 and (n) = 0.5 in Fig. 10 for which we obtain 
a g-wave, and higher order s-wave, respectively. In the 
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